Ultra-long distance interaction between spin qubits 
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We describe a method for implementing deterministic quantum gates between two spin qubits 
separated by centimeters. Qubits defined by the singlet and triplet states of two exchange coupled 
quantum dots have recently been shown to possess long coherence times. When the effective nuclear 
fields in the two asymmetric quantum dots are different, total spin will no longer be a good quantum 
number and there will be a large electric dipole coupling between the two qubit states: we show that 
when such a double quantum dot qubit is embedded in a superconducting microstrip cavity, then 
the strong coupling regime of cavity quantum electrodynamics lies within reach. Virtual photons in 
a common cavity mode could mediate coherent interactions between two distant qubits embedded 
in the same structure; the range of this two-qubit interaction is determined by the wavelength of 
the microwave transition. 



Introduction. Experimental realization of conditional 
quantum dynamics of two isolated solid-state quantum 
systems has become a holy grail of mesoscopic physics re- 
search, due to its potential implications for scalable quan- 
tum information processing. The majority of theoretical 
proposals aimed at this goal is based on nearest-neighbor 
interactions such as the Heisenberg exchange coupling be- 
tween quantum dot (QD) spins 1]. However, to achieve 
lower accuracy thresholds for quantum error correction, 
the implementation of coherent long-range interactions 
between two qubits is highly desirable [|J • Optical dipolc- 
dipole interactions 0] , capacitive coupling |i( and optical 
cavity-mediated interactions between spins could be 
used to realize controlled quantum gate operations on 
length-scales comparable to optical wavelengths: these 
mechanisms may then enable coherent interactions be- 
tween a limited number (< 10) of QD spins. 

In this Rapid Communication, we show how ultra- 
long range coherent interactions between two spin qubits 
separated by centimeters can be mediated. Our pro- 
posal is motivated by two recent remarkable experimen- 
tal achievements: 1) realization of circuit-QED using a 
Josephson-junction charge qubit strongly coupled to a su- 
perconducting (SC) microstrip cavity |5|; and 2) demon- 
stration that the singlet-triplet subspace of a double-QD 
structure constitutes a promising qubit exhibiting coher- 
ence times exceeding 10 /is @, 0|- We show here that due 
to the presence of magnetic field gradients caused by par- 
tially polarized QD nuclear spin ensembles, it is possible 
to induce a large electric-dipole coupling between singlet 
(S) and triplet (To) states by adjusting an external gate 
voltage 8] . Since the energy of the S-Xb transition is in 
the microwave range, it is possible to use SC microstrip 
cavities with a length (L) equal to the transition wave- 
length (A) and a cavity-volume V cav ~ 10 _8 A 3 to me- 
diate interactions between two qubits embedded in the 
same cavity via virtual microwave photon exchange. A 
distinguishing feature of our proposal is the large separa- 
tion between the length scales determining single qubit 
control, determined either by fabrication (~ 50 nm) or 
optical wavelength (~ 1/im), and that of two-qubit in- 




FIG. 1: (Color online) The proposed setup with two vertically 
coupled double QDs next to a superconducting micro-strip 
cavity. The topmost QD serves as a marker. 



teractions, ultimately determined by the wavelength cor- 
responding to the (adjustable) S-Tq qubit transition. As 
in Ref. 0, the qubits here are coupled via their electric 
dipole moment: however, the use of a cavity does away 
with the 1 /r 3 decay of the dipolar interaction that limits 
the range of the coupling in Ref. 0. In contrast to the 
scheme introduced in Ref. [fj which couples spins via their 
magnetic dipole moment, our proposal does not require 
the use of electron spin resonance (ESR). 

Fig. [T] shows the structure that we envision: the micro- 
strip cavity is defined by a wavelength- long center SC 
strip separated from the ground planes by ~ 100 nm. 
The whole structure is deposited on a Molecular Beam 
Epitaxy grown GaAs wafer containing a stack of two self- 
assembled QDs that are tunnel coupled and buried ~ 
100 nm below the surface. Finally, ~ 50 nm below the 
lowest QD layer lies either an n-doped 20 nm GaAs layer 
or a modulation-doped quantum well. The ohmic contact 
to this bottom electron reservoir allows for applying a 
gate voltage Vgato that is used to inject single electrons 
into each QD dctcrministically and to bring the electronic 
states of the two QDs in and out of resonance [1(3] ■ Even 
though the self-assembled QDs in the first layer nucleate 
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at random locations, the QDs of the second layer have 
a very high likelihood for nucleating directly above the 
QDs of the first layer. It has been shown that Atomic 
Force Microscopy can be used to determine the position 
of stacked QDs with a spatial resolution of 25 nm . In 
order to apply independent gate voltages to two double 
QDs embedded in the same cavity, it will be necessary 
to wafer-fuse two separate samples before depositing the 
SC thin layers 0. 

Coupling mechanism. Coupling the spin singlet |5) 
and triplet \Tq) in a double QD (the two qubit states) via 
the emission or absorption of a cavity photon requires a 
sufficiently strong electric dipole transition between the 
two states [l3[. The key question is under what condi- 
tions it is possible to obtain an electric dipole transition 
in a double QD. First, we remark that the two QDs need 
to be coupled via inter-dot tunneling (with a tunneling 
energy t) for a non-zero dipole matrix element; indeed, 
we find below in Eq. ([2]) that the matrix element is pro- 
portional to the singlet-triplet energy splitting, i.e. the 
exchange energy Joe t 2 QJ. 

Provided that \t\ > 0, there are still two indepen- 
dent symmetries that can prevent electrical dipole tran- 
sitions. One of the two symmetries derives from the 
spin-conserving nature of the electron-photon interac- 
tion. The spin singlet |5) and triplets \T ,±) are eigen- 
states of the total spin with different spin quantum num- 
bers 5 = and 5 = 1 and cannot be transformed into 
each other by the emission or absorption of a photon 
which changes only the orbital angular momentum. How- 
ever, |5) and \Tq) are mixed, and thus the spin selection 
rules broken, by the presence of a magnetic field that 
is inhomogeneous on the scale of the inter-dot distance. 
Such a field inhomogeneity Sh is usually unavoidable in 
the form the Overhauser field due to the hyperfine cou- 
pling of the electron spin to the surrounding nuclear spins 
in the QD material. The dipole matrix element given be- 
low in Eq. @ is indeed proportional to Sh. 

The second problem to be overcome if dipole transi- 
tions are to occur between |5) and |T ) in a double QD 
is the orbital symmetry that exchanges the two QDs. 
The effect of tunneling on 15) consists in the admix- 
ture of the \S(1, 1)) singlet (one electron in each QD) 
with the states 15(2,0)) and 15(0,2)) that involve two 
electrons on the same QD; for symmetric QDs, this 
admixture is restricted to the symmetric combination 
\D+) = (|5(2,0)) + |5(0,2)))/V2 of doubly occupied 
states on the two QDs, while the electric dipole Hamil- 
tonian has odd parity (it can be represented in terms 
of the momentum or the position operator, both having 
odd parity) and thus couples the singlet (even parity) 
exclusively to the antisymmetric combination |_D_ ) = 
(|5(2,0)) - |5(0,2)))/V^, and thus not to the "qubit" 
singlet. The mirror symmetry between the QDs is bro- 
ken if the QDs are electrically biased, thus detuning their 
single-electron levels by an energy e. Roughly speaking, 
the electric bias e creates a situation with a mobile charge 
that allows for an electric dipole moment which is absent 
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FIG. 2: (Color online) Energy of the two-electron states in a 
double QD as a function of the inter-dot detuning e around 
the resonance e = U, indicated by a dotted vertical line. The 
parameters chosen for the plot are U = lOmeV, t — 0.1 meV, 
Sh — 0.15 meV, and p/is-B = 1 meV. 



in the unbiased double QD (e = 0). We plot the two- 
electron spectrum in a double QD as a function of s in 
Fig. [2 We expect that the dipole matrix element is pro- 
portional to e. 

Having discussed the underlying physical considera- 
tions, we proceed with the key results of the paper and 
defer the derivation of qubit-cavity coupling strength g 
to the last part of the Letter. In the presence of electric 
(or magnetic) dipole coupling between the states |5) and 
\Tq) that define our qubit, the Hamiltonian is 



H = -a z + ga x (a + a f ) , 



(1) 



with e the 5- To splitting and, for electric dipole coupling, 



g = eaEo 



J 



e {Sh/2) 



Tilj U 2 -e 2 - (Sh/2) 2 



(2) 



The vacuum value of the electric field is given by Eq = 
Kms/^ = \/huj/2e a eLd 2 [III, where L is the length of 
the center SC, d is its separation from the ground SC 
planes, and e is the effective dielectric constant seen by 
the cavity mode. The operators a< and a create and 
annihilate a cavity photon with frequency w/2ir. The 
energy denominator in Eq. @ arises from the admixture 
of the 5(1, 1) and Tq states with the doubly occupied 
states 5(2, 0) and 5(0, 2) that are separated in energy by 
U ± e ± Sh/2. The details of the derivation of Eqs. Q 
and ^ will be given further below. 

Two- qubit coupling. We now turn to the situation of 
two double QDs coupled to the same cavity, as shown 
in Fig. [T] By introducing the rotating wave approxima- 
tion in Eq. ([T]) and eliminating the cavity mode using a 
Schrieffer- Wolff (SW) transformation 0,13], we obtain 



8=1,2 



eff I cr+ o_ + a_ a + J , (6) 



with the effective qubit-qubit coupling parameter g e g = 
<7i52 [l/(^i — + 1/(^2 — Jt-uj)], and the Stark-shifted 
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single-qubit splitting, ii/2 = e,/2 + g 2 ((n) + 1/2) /(e t - 
fuS), where (n) = (a^a) denotes the number of photons 
in the cavity. 

Qubit- cavity coupling. To derive Eqs. ([!} and ([5]), we 
consider the Hamiltonian of a single qubit in the cavity, 



spanned by |T ), \S), \D+), and \D. 



H 




Sh/2 

2f 

2t U e 

e U 



(6) 



H — H c \ + H cav + i/dip , 



(4) 



describing the electronic degrees of freedom, the cavity 
field £f C av = hu(a^a + 1/2), and the electric dipole cou- 
pling between the qubit and the cavity. In the first step 
of our derivation, we concentrate on the electronic part: 
we write H e \ = Hp + Hj- + H int where 



a = L,R 



(5) 



denotes the 
lowest-energy 

Ht = 

for electron 

flint = UJ2 a 



single-electron Hamiltonian of the 
orbital on each separate QD, 



■r.[ 



accounts 

tunneling between the QDs, and 

l R c ai Ca l c ai Ca l describes the Coulomb 
interaction between two electrons occupying the same 
QD. In Eq. (O, e — El — £r is the asymmetry of 
the double QD and B LtR = B ± SB/ 2 the effective 
magnetic field for an electron on the left (right) QD. 
The presence of nuclear spins in the QDs gives rise to 
inhomogeneous effective magnetic (Overhauser) fields 
Bl.r — J2iAiI LR , where I l L R is the ith nuclear spin 
in the QD L, R, and Ai is the corresponding hypcrfinc 
coupling constant 17(. The operators cj^ (c aa ) create 
(annihilate) an electron in an orthogonalizcd Wannicr 

orbital $ LjR = (ifiLM. ~ 7'Ar,l)/v /i _ 25 7 + 7 2 in tnc 
QD a = L,R with spin a =|,J,, where S = ((^lIv'-r) 
denotes the overlap integral between the left and right 
unnormalized orbitals and 7 = (1 — y/l — S 2 )/S. In 
parabolic QDs, the wavefunctions ip a are Gaussian. 

A low energy two-electron double QD where only 
the ground orbital state on each QD can be occu- 
pied, has six possible states: the three spin triplets 



TO 



L1 C Rl 



C Ll c R1 



10), \T a 



-■Ra 



|0) 



i°~ =T>J.) with S z 



glets \S) = |-S(1,1)> = ^( 



0, ±1, and the three spin sin- 
|0), and 



LI Rl 



\D±) = ^ (c[^ Li ± 4 T 4J |0), all with S = S z = 0. 

Here, D± are linear combinations of the states with dou- 
ble occupation of a QD and |0) is the state with no elec- 
trons. 

We choose a coordinate system such that the z axis is 
along the homogeneous part of the field B and decompose 
the difference field into its longitudinal and transverse 
parts, <5B = <5B 2 + 5Bj_. We assume 5B± <c B z which 
ensures that the spin-polarized states \T a ) are decoupled 
from the remaining four states [l8l |. We can then write 
the Hamiltonian as the four-by-four matrix in the basis 



where we have introduced the relative Zeeman energy 
8h = g{iB8B z between the dots. For e = 0, the \D-) 
state completely decouples because its orbital symmetry 
forbids any coupling to the other singlets, while a cou- 
pling to the triplet is impossible due to spin conservation. 

In the weak tunneling regime t -C U — e we can 
eliminate \D±) by means of a SW transformation [l6[ 
H = er s He s ~ H + [H T , S] /2, with S = and 
flo = flo + flint- The terms of order Ht are cancelled 
in A because we choose S such that [Hq, S] = —Ht, 



S 



At 



{U 2 



2Sh 2 (U 2 +e 2 ) 



s 
-s 



(7) 



with 



_ / Sh(U 2 +e 2 - (5h/2) 2 ) -2ShUe 
S ~ \ 2U{U 2 -e 2 - {Sh/2) 2 ) -2e(U 2 - e 2 + [Sh/2) 2 ) 

We assume here that we are in the regime U 3> t, Sh and 
allow e to lie in the whole range < e < U. The effect 
of the SW transformation is to separate the states with 
single occupation from \D±) within the lowest order in 
t/(U — e) in the Hamiltonian, 



H 



H s 
H D 



with the exchange coupling 



At 2 U(U 2 



(Sh/2) 2 ) 



4. ( 10 ) 



(U 2 - e 2 ) 2 - 2(5h/2) 2 (U 2 + e 2 ) + (Sh/2) 

and the effective relative field 

i Wl, J(U 2 +e 2 -(5h/2) 2 ) , 
6h = 5h 1 1 - WW -e 2 - (Sh/2) 2 ) ) ' 



Note that for e = Sh = 0, Eq. (fTO]) reduces to the famil- 
iar expression J = At 2 /U . From Eq. (|9|), we obtain the 
singlet and triplet eigenenergies 



£± = i ( ■-./■; ! A/ / - 



2 2 



(12) 



Dipole matrix element. Optical transitions conserve 
spin, therefore transitions between the singlets \S), \D±) 
and the triplets |Tj) (i = 0, |,4) are forbidden. How- 
ever, the presence of an inhomogeneous magnetic field 
SB breaks spin symmetry and allows for electric dipole 
transitions. The dipole coupling to a single cavity mode 
is described by [lj| 



fld 



ip 



6 A e 

— A -p = . 

m m \2eoeVuj 



1/2 



e ■ p (a + a' 



(13) 
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where V = Ld 2 , Aq denotes the vector potential at r = 
and (a) creates (annihilates) a cavity photon with 
frequency u>, described by fl c „. In the following, we 
determine the dipolc matrix element 

■9=--f^77^) 1/2 ^l e -P^)' ( 14 ) 
m \ZSqsVlu J 

where \S) and \Tq) are the eigenstates of Hs in Eq. ([9]). 

In the single-particle eigenbasis of Hp + Hx, 
i.e., the bonding and antibonding orbitals 

|*±,*) = (cL ± cL)|0)/V2, and using (n\p\m) = 
-im{n\[H D + H T ,x\\m}/h, we find ($±|p|$±) = 
and ($_|p|$+) = (-imt/ h)[{tp L \x\tp L ) - (<p R \x\ip R ) + 
2ilm((pi J \x.\ipfi}]/2\/l — S 2 . The last term can have a 
nonzero component perpendicular to x due to orbital 
diamagnetism, but it turns out that this real contri- 
bution to the single-particle matrix element does not 
contribute to g. We find |p 2; | < i>+) = imta/h, while 
the other components have vanishing imaginary parts. 
The electric dipole moment of the double QD is thus 
directed vertically in Fig. [1] and couples to the vertical 
component of the cavity field which is increased by 
positioning the QDs close to the microstrip. 

The only nonvanishing two-electron matrix clement be- 
tween the unperturbed states \S), \Tq), and \D±) is 

(D-\p x \S) = 2ilm(<f>-\p x \$ + ) = 2imta/Ti. (15) 

We can transform -f/dip into the new basis using the same 
SW transformation -f/dip — ^dip + [-ffdip, S]. In the sub- 
space spanned by the transformed states \S) and |To), we 
obtain for the momentum operator in i?dip 

. _ iam eJSh/2 ( -1 \ 

Px ~ h U 2 -e 2 - (5h/2) 2 \l ) ■ {lb > 

Transforming Eq. (TIT)]) into the eigenbasis of Eq. (0) is a 
rotation in the (pseudo)-xz plane thus leaving the dipole 
Hamiltonian Eq. ([16)) , having the form of a pseudo-field 
in y-direction, invariant up to a phase factor. With 
Eq. p4p , the qubit-cavity Hamiltonian Eq. ((3]) in the log- 
ical subspace of the singlet-triplet qubit takes the form 
of Eq. (fTJ) with the coupling constant Eq. ([2]). Close 



to resonance, we can replace Tiuo by e = e+ — S- = 
VJ 2 + Sh 2 in Eq. ©. For htu ps 0.1 meV, e ~ 13 
(GaAs), V = 1 cm(100nm) 2 , we arrive at a vacuum field 
of Eq ps 25V/m. For self-assembled QDs, we estimate 
the dot-distance to be of the order of a « 10 nm, there- 
fore eaE /h ~ 0.25 ^eV ~ 65 MHz. 

For further discussion, we introduce A 2 = U 2 — e 2 
and work in the regime of a strongly biased QD pair, 
A <C U. We can envision a hierarchy of energy scales 
Sh, J <C A <C U, e.g., J w Sh w 0.1 meV (for recent ex- 
perimental results on optically generated nuclear polar- 
ization in QDs, see Ref.Q), A > 1 meV, and U ~ 10 meV. 
In this case Sh Sh, and near resonance g « eaE$. 
Thus, g ^> u>/Q, 7, provided that the cavity quality factor 
Q > 10 4 and the spin decoherence rate 7 < 10 7 s _1 . The 
SW transformation can be applied if t/V A 2 — Sh 2 <C 1: 
hence the "resonance" A — > Sh where formally g — ► eai?o 
is not within the regime of the validity of our result. 
If J <C 5/i <C A -c U, we can simultaneously sat- 
isfy USh/A w 1 (3 fa ea£: ) and i/ A < L In tlic 
regime e, <5/i <C J7, we obtain Sh w 5/i and therefore 
5 ~ eaE s Sh/U 2 -C ea_E - 

Conclusions. A cavity-doublc-QD coupling strength 
of g ~ 65 MHz implies that it is possible to implement 
two-qubit gates on timescales ~ 10 ns: while this is al- 
ready three orders of magnitude shorter than the spin 
(memory) decoherence time, an important open question 
that needs to be addressed is the gate errors. The fact 
that total spin is not a good quantum number during the 
gate operation when the dipole-coupling of \S) and \Tq) 
states is on, most likely introduces additional phonon or 
charge fluctuation mediated decoherence. However, even 
in the case of relatively strong charge decoherence, one 
could still envision using the cavity-mediated coupling as 
a source of entanglement that can be distilled and used 
for remote gate operations. We also emphasize that the 
observation of conditional quantum dynamics between 
two spins with a macroscopic separation would itself be 
an exciting goal for the emerging field of spintronics. 
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